8689 


' 


✓ 


ORC  68-23 
SEPTEMBER  1968 


SOME  RESULTS  FOR 

INFINITE  SERVER  POISSON  QUEUES 


SHEiOON  M.  ROSS  AND  MARK  BROWN 


OPERATIONS 

RESEARCH 

rfeNTER 


A 


ThS  docuznan*  bca 
for  public 

disMbqRoo  It 


D  D  C 

ropor-'T  ^fnl 
i  NOV,  4  1968 

JkiaKirmdli) 

f  8. 


Reproduced  by  the? 

CLEARINGHOUSE 

lor  federal  Scientific  &  Technical 
Information  Springfield  Va  22151 


UNIVERSITY  OF  CALIFORNIA  •  BERKELEY 


AD  6 


SOME  RESULTS  FOR  INFINITE  SERVER  POISSON  QUEUES 


by 


Sheldon  M.  Ross 

Department  of  Industrial  Engineering 
and  Operations  Research 
University  of  California,  Berkeley 

and 

Mark  Brown 

Department  of  Operations  Research 
Cornell  University 
Ithaca,  New  York 


September  1968  ORC  68-23 

This  research  has  been  partially  supported  by  the  U.S.  Army  Research 
Of f ice-Durham  under  Contract  DA-31-124-ARO-D-331  and  the  National 
Science  Foundation  under  Grant  GP-8695  with  the  University  of 
California.  Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government. 


ABSTRACT 


A  generalization  of  the  M/G/®  queueing  system  with 

batch  arrivals  to  one  with  time  dependent  arrival  rates, 

service  times,  and  batch  size  distributions  is  considered. 

It  is  shown  that  both  W(t)  ,  the  number  of  people  being 

served  at  t  ,  and  S(t)  ,  the  number  of  people  who  have 

completed  service  by  t  ,  are  distributed  as  compound  Poisson 

laws.  The  distributions  of  the  traffic  time  average 
-1  T 

T  f  W(t)dt  and  the  occupation  time  0(t)  (the  amount  of 
0 

tine  past  t  until  the  system  becomes  empty,  under  the 
assumption  that  no  new  customers  are  served  after  t)  are 
also  derived. 

The  limiting  proportion  of  busy  time  and  the  asymptotic 
behavior  of  the  traffic  time  average  are  also  discussed  in 
the  time  homogeneous  case. 
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0.  Introduction  and  Summary 

We  consider  a  queueing  model  in  which  arrivals  occur  according  to  a 
nonhomogeneous  Poisson  Process  in  batches  of  varying  size,  and  in  which  a 
customer  is  served  immediately  upon  arrival  by  one  of  an  infinite  number  of 
servers. 

We  allow  for  the  possibility  that  both  the  batch  size  and  service  time 
distributions  might  depend  on  the  arrival  time  and  thus  denote  by  P^fr)  ,  the 
probability  that  a  batch  arriving  at  time  t  will  contain  r  customers 
(r  >  1)  ,  and  by  Gt()  ,  the  service  time  distribution  of  a  customer  arriving 
at  time  t  .  We  also  let  m(t)  denote  the  mean  value  function  of  the  Poisson 
Process  of  arrivals.  This  system  is  thus  the  generalization  of  M/G/°°  with 
batch  arrivals  to  time  dependent  arrival  rates,  service  times  and  bat:h  size 
distributions. 

In  the  first  section  we  show  that  both  W(t)  ,  the  number  of  customers  being 
served  at  time  t  ,  and  S(t)  ,  the  number  of  customers  who  have  completed 
service  by  time  t  ,  are  distributed  as  compound  Poisson  laws.  In  the  second 
section  we  derive,  in  the  time  homogeneous  case  (Gt()  ■  G()  ,  Pt()  ■  P()  , 
m(t)  ■  At)  ,  the  limiting  proportion  of  time  that  the  system  is  nonempty. 

In  the  third  section  we  derive  the  distribution  of  the  occupation  time 
0 ( t )  ;  where  0(t)  is  defined  as  the  amount  of  time  past  t  until  the  system 
becomes  empty,  under  the  assumption  that  no  new  customers  are  served  after  time  t 

In  the  fourth  section  we  derive  the  distribution  of  the  traffic  time  average 
-1  T 

T  f  W(t)dt  ,  and  its  asymptotic  behavior  is  discussed  in  the  time  homogeneous 

0 


case. 


2 


In  [7]  Shanbag  considered  a  special  case  of  the  above  model,  G  ■  G  , 

?tO  ■  P()  ,  and  by  solving  a  differential  equation,  derived  the  joint  generating 
function  of  W(t)  and  S(t)  .  His  method,  however,  does  not  seem  applicable  to 
the  present  model  unless  some  conditions  (such  as  t-continuity)  are  placed  on 
Pt()  and  Gt()  . 

Benes  (rf.  [6],  p.  123)  has  previously  obtained  the  distribution  of  the 
traffic  time  average  for  the  case  M/M/®  ,  and  in  a  more  recent  paper  [5]  Rao 
has  generalized  this  to  the  case  CG.I/G/®  where  CG.I  stands  for  any  stationary 
stream  of  random  jumps  (batch  arrivals)  for  which  the  times  between  successive 
jumps  are  independent  and  identically  distributed.  His  results  thus  include 
M/G/®  with  batch  arrivals  as  a  special  case.  The  method  employed  in  the  present 
paper  differs  from  those  used  in  the  above  papers. 


« 


t 
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1.  Distribution  of  W(t)  and  S(t) 

Throughout  this  paper  we  shall  assume  Px(r)  is  a  measurable  function  of 
x  for  all  r  end  that  G^Ct  -  x)  is  a  measurable  function  of  x  for  all  t  . 
We  shall  suppose  that  the  process  first  begins  at  t  ■  0  ,  and  we  let 

B(t) 

B(t)  be  the  number  of  batches  which  have  arrived  by  time  t  .  Then  W(t)  *  £ 

i-1 

where  y^  denotes  the  number  of  arrivals  from  the  iC^  batch  that  are  being 
served  at  time  t  .  Thus 

P{W(t)  -  k}  -  l  e"m(t)  pj[  yt  -  k  |  B(c)  -  n|  .  (1) 


Now  conditional  on  B(t)  *  n  ,  the  (unordered)  arrival  times  of  the  batches 

are  distributed  as  a  sample  of  independent  and  Identically  distributed  (i.i.d.) 

m(x)/m(t)  x  <  t 


x  >  t 


random  variables  with  a  distribution  given  by  F(x)  * 

B(t) 

Thus  conditional  on  B(t)  «  n  ,  £  y.  is  distributed  as  the  sum  of  n  J.i.d. 

1  1 

random  variables  Z^,  each  having  probability  distribution 


P(2l  '  J1  ’  {  X  ?x(,>(j)w  -  G*(t  -  ^  <Gx(t 


-  x))r  ^dm(x)  (2) 


Thus, 


P{W(t)  -  k}  -  l  e"m(t)  P{ZX  +  ...  +  Zn  -  k}  (3) 

n 


and  so  W(t)  has  a  compound  Poisson  distribution  with  Poisson  parameter  m(t) 

B(t) 

and  with  jumps  distributed  according  to  (2)  —  i.e.,  W(t)  -  l  where 

are  i.i.d.  according  to  (2)  and  are  independent  of  B(t)  .  The  probability 


- JdC 


generating  function  *w(t)(s)  =  E(sW(t))  ls  given  by 


Vt)(s)  -  «*p}jI0  <sJ  -  1>  {lmi  Px(r)(j)d  -  Gx(t  -  x))J(Gx(t  -  x))r~^dm(x)  j  . 

A  similar  analysis  may  be  done  to  show  that  S(t)  has  a  compound  Poisson 
distribution  with  Poisson  parameter  m(t)  and  Jump  distribution  V  where 

P  V  ^  "  “o* md)  r^j  ~  d  “  V*  "  x^r  ^dm(x)  ,  j  >  0  .  (5) 


'i.'Srw*,... 
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2.  Limiting  Proportion  of  Busy  Time  (Homogeneous  Case) 

In  this  section  we  suppose  that  m(t)  "  At  ,  G  ■  G  ,  and  P  ■  P  . 

From  (2)  and  (3)  it  follows  that 

P{W(t)  *  0}  ■  exp  | -A  f  \  (1  -  Gr(x))  P(r)dx|  ->  e  as  t  -*■  00  ,  where  (6) 

00 

M  mf  l  (1  -  Gr(x))  P(r)dx  .+ 

0  r 

Now  as  time  passes,  there  will  be  periods  of  time  at  which  the  queue  is  empty 
which  wixl  alternate  with  periods  at  which  the  queue  is  busy.  Let  A^  be  the 
length  of  the  i**1  empty  period  and  the  length  of  the  iC^  busy  periods. 

OD  00 

The  sequences  {A.}  and  (D, }  are  independent  renewal  sequences,  and  thus 
I  1  1  l 

the  sequence  fAi»Di*A2,D2’  ...}  is  an  alternating  renewal  sequence.  Then  it 
is  known  (see  [4])  that  P{Queue  is  empty  at  t}  -*•  -  ~jy(+ED  as  t  00  , 

and  thus  from  (6)  we  have  that 


ED 


tt 


(7) 


Let  +  Di  »  ch«n  (C^}  is  a  renewal  sequence.  Let  N(t)  denote  the 

number  of  C-renevals  up  to  time  t  .  Since 

N(t)  a. s .  _1  -AM 
t  -*■  EC  * 

*If  r  denotes  a  random  batch  size,  and  Y^,  Yf  the  service  times  then 

M  ■  E[Max(Y. ,  ....  Y  )]  and  this  is  finite  if  f  (1  -  G(y)dy)  and  £rP(r)  are 
1  r  0 

finite. 

++Equation  (7)  vaa  derived  by  a  different  method  in  [7]. 


6 


and 


t  (see  [2])  it  follows  from  the 


strong  law  of  large  numbers  that 


,  N(t) 

\  l  At  *i 
1  1 


a.s.  -AM 

e  as  t  -►  *  , 


(8) 


and  from  Wald's  Fundamental  Equation  of  Sequential  Analysis  th*t 


r  "(o  1 

:[1/c  |  AiJ 


-AM 

e  f  j  t  -*■ 


(9) 


Now  let  A(t)  ”  amount  of  time  the  queue  is  empty  up  to  time  t 
D(t)  ■  amount  of  time  the  queue  is  busy  up  to  time  t  . 

N(t) 

Now  |  A(t)  -  i  At  |  <  AN(t)+1  ,  and  since  1/t  E[AN(t)+1]  -  0  and 
*11(0+1  0  with  Probl«®  1  as  t  -*•  •  we  have  by  (8)  and  (9)  that 


Aft)  a.s.  -AM 
j  e  as  t 


and 


rfA(t)1  -AM 

EL  t  J  e  *8  t  *  • 


Also  since  D(t)  •  t  -  L(t)  we  have  that 


(10) 


and 


D(t)  a^s.  1  e-AM 


*  1  - 


(ID 


V 


M, 

7 


It  can  also  be  shown  that  A(t)  and  D(t)  (suitably  normalized)  both 
»  having  limiting  normal  distribution  (see  [8]). 


ft 


* 


* 


The  occupation  time  0(t)  is  defined  as  the  amount  of  time  past  t  until 


the  system  becomes  empty  when  no  new  customers  are  served  after  time  t  . 


We  say  that  a  batch  is  served  when  all  members  of  that  batch  have  been 


served.  Now  the  time  points  at  which  batches  being  served  at  time  t  arrived 
may  be  shown  (see  [3],p.  497  or  [l],p.4)  to  form  a  nonhomogeneous  Poisson  Process 

y 

with  mean  value  function  m(y)  -  f  (1  -  G  (t  -  x)dm(x)  ,  y  <  t  ,  where 

0  X 

G  (a)  ■  \  P  (r)Gr(a)  ;  and  thus  given  that  there  are  n  batches  being  served 
r 


at  t  their  (unordered)  arrival  times  have  the  same  distribution  as  an  i.i.d. 


sample  from 


y  t 

F(y)  -  f  (1  -  G  (t  -  x)dm(x) /  J  (1  -  G  (t  -  x)dm(x)  y  <  t 
0  X  0  X 


and  so 


0t(x)  =  P(0(t)  <  x) 


r  /  'A(t  +  x  -  y)  -  5  -  y>  \n 

-  I  P{W(t)  -  nH/-* - - - 1 -  dF(y) 

n  \  0  1  -  Gy(t  -  y)  /  (12) 


-  ♦y(t)(  •/*  (Gy(t  +  x  -  y)  -  Gy(t  -  y))  dm(y)/  /  (1  -  Gy(t  -  y)dm(y)j  . 


9 


4.  Traffic  Time  Average 

1  T 

In  order  to  obtain  the  distribution  of  W  -  —  f  W(t)dt  we  first  note 

T  0 


/  W(t)dt 
0 


B(T)  ri 

[  l  Min  (x  ,  T  -  t  ) 
i-1  j=l  1J  1 


where 


*  arrival  time  of  i  batch 


*  number  in  i  batch 


•  service  time  of  member  of  i**1  batch 


and  thus 


T  B(T) 

f  W(t)dt  -  l  L 
0  i-1  1 


T  £tl 

where  Is  the  sum  at  T  of  all  the  service  times  of  members  of  the  i 

batch.  It  thus  follows  as  in  Section  1  that 


T  B(T) 

/  W(t)dt  -  I  R.  , 

A  J  l  ^ 


where  R^  are  i.i.d.  independent  of  B(T) 


and  where 


where 


p(Ri <  *>  ■  {  im  l  px(r)Gx^ 
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ir 


a  <  T  -  x  , 


and 


(a) 


is  the  r-fold  convolution. 


Letting 


f  e^UadG  (a)  ,  we  have  that 
0  x,i 


Vu)  5 

T 


IT 


(17) 
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5.  Homogeneous  Case 

We  suppose  that  Gf  *  C  ,  c  P  ,  and  M(t)  =  At  ;  also  pg  = J rP(r)  , 

2  2 
p  _  *  Vr  P(r)  ,  p  ■  f  xdG(x)  ,  and  p  =■-  f  x  dG(x)  are  .ill  assumed  finite. 

B2  L  G  G2 


Let  be  the  sum  of  the  service  times  of  members  of  the  i 


th 


batch,  and 


1Ct  WL  "  UB  UG  *  3nd  W  2 

L 


“b  +  V  “c 


•  |o2  =  y  (x  -  UG)2dG(x)j  .  Let 


_  /  B(T)  -  \ 

sx  ■  a  (1/T  I  Li  -  x"l) 

ST  "  *  (1/T  [  Li  -  X^l)  • 

Now,  Var  (s*  -  STJ  -  a(pl  -  EL*J  +  A^o2  -  Var  L*  J  -  0  as  T  -  •  .  Also 
e|s^,  -  S^j  0  .  Thus,  converging  in  distribution  implies  that  ST  also 

converges  in  distribution  to  the  same  limit  law.  Now 


*  *(t)  «=  exp  j  AT  ($L(t//r)  -  1)  -  it  /T  pl 


and 


♦L(t/^)  -  1  +  it  ^//T  -  t*  p  2/2T  +  oCT"1)  , 

L 


Implying  that 


as  T 


(19) 


Thus, 
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ft  (WT  -  ^PL)  l  Normal  (0  ,  Au  2)  .  (20) 

Now,  let  N(a,b)  be  the  number  of  customers  arriving  in  (a,b)  ,  and  let 

N(t)  -  N(0, t)  . 

« 

Lemma  1 : 

W(t)/N(t)  a.'>S'  0  as  t  -*■  00  . 

Proof : 

W(t)/N(t)  <  N (t  -  n,t)/N(t)  +  W(t ,n)/N(t  -  n)  ,  where  W(t,n)  denotes  the 
number  of  customers  arriving  in  (0,t  -  n]  whose  service  time  is  greater  than 
n  .  Thus, 


lim  W(t)/N(t)  <  1  -  G(n)  a.s.  for  all  n  . 

t-**> 


Q. E.D. 


Theorem  1: 


T 

1/T  f  W(t)dt  V  Xy.  as  T  -*■  ®  . 
0  L 


Proof : 

Suppose  first  that  service  times  are  bounded,  i.e.,  G(M)  ■  1  for  some 
M  <  *  ,  and  let  {X^  ,  i  ■  1,  ....  W(T)}  be  the  service  times  of  customers  being 
served  at  T  .  Then 

B(T)  B(T)  _  W(T) 

1/N(T)  l  Lj  -  1/N(T)  l  L*  <  J  Xj  *  0  by 


4 
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Lemma  1.  Thus, 


B  (T) 


1/T  /  »(t)dc 


and  so  the  result  follows  In  the  bounded  case. 

Now  suppose  only  that  pD  and  u  are  finite.  Let  W^ft)  denote  the 

B  G 

number  of  customers  being  served  at  t  whose  ser  ice  time  at  t  is  less  than 
M  .  A1  so,  let  N(t)  denote  the  number  of  batches  arriving  in  (0,t)  having 
a  member  whose  service  time  is  greater  than  or  equal  to  M  ,  and  let  be 

the  sum  of  the  service  times  of  the  i^  such  batch.  Then 


T  T  N(T) 

1/T  f  W(t)dt  -  1/T  f  W  (t)d t  <  1/T  "  L 

0  0  ‘  i  1 


(21) 


AEL  (l  -  l  P(r)Gr(M)] 
-*=°  1 V  r  / 


a 
T-~ 


(22) 


where  the  convergence  follows  from  the  fact  that  N(t)  is  a  Poisson  Process  with 


mean 


value  function  \t^l  -  £  P(r)Gr(M)j  .  Now 


EL±  -  l  P(r )/(y:  +  ...  +  yr)dG(yi)  ...  dG(yr)/l  -  J  P(r)Gr(M) 


(23) 


r  * 

M 


where 


M  ■  (Max  (y^,  ...,  y  )  >  M). 


Thus,  combining  (22)  and  (23)  we  arrive  at 


1/T  f  (W(t)  -  w”(t))dt  <  \l  P(r)  f  (y  ,  ....  y  )dG(y  ) 

0  r  * 

M 


.  dG(yf)  a.s. 


(24) 


as  T  -►  ®  .  However,  the  right-hand  side  of  (24)  goes  to  zero  as  M  -►  °° 
#  Pg  and  p^  are  finite),  and  so  the  result  follows  from  the  bounded  case 

either  pg  and  is  infinite,  the  result  follows  from  truncation. 


(since 

If 
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